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Selecting specific initial configuration using spectator neutrons in U+U collisions 
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We present a method to select events with specific initial configuration, namely body-tip, in heavy- 
ion collisions using deformed Uranium nuclei. We propose to use asymmetry in spectator neutron 
numbers to filter out these body-tip events from the unbiased configurations in U+U collisions. 
We have used a variable S v to differentiate between the body-tip and unbiased configurations. We 
have calculated the 2 nd order azimuthal anisotropy, namely elliptic flow ( 02 ), for this body-tip 
configuration in the framework of a transport model and found it to be consistently lower compared 
to that in unbiased configurations as we expected. The purity of selecting such events in a real 
experiment is also discussed. 


I. INTRODUCTION 

In a central heavy-ion collisions with spherical nuclei such 
as Au or Pb, the initial overlap region is always circular. 
In the U+U collisions, however, the initial overlap re¬ 
gion can acquire different configurations owing to the de¬ 
formed shape of the Uranium nuclei [l;j. In Ref. [lj, the el¬ 
liptic flow (v 2 ) has been studied and found to be strongly 
correlated with the different configurations of the initial 
overlap region. Moreover, in heavy-ion collisions, the en¬ 
ergetic spectator protons can produce a strong magnetic 
field reaching eB y 0. Such a strong magnetic field 

can give rise to chiral mag netic effect (CME) and chiral 
separation effect (CSE) [JJ3|. However, to observe these 
phenomena in real data, the azimuthal anisotropy (+ 2 ) 
has to be minimised as it acts as a background 3 to these 
processes. Therefore, the events which have very high 
magnetic field and low azimuthal anisotropy (v 2 ) are the 
perfect candidates for the CME. In a central Au+Au or 
Pb+Pb collisions, the azimuthal anisotropy is very low, 
but the magnetic field is also low due to less number of 
spectator protons. In a non-central Au+Au or Pb+Pb 
collisions, although the magnetic field is comparatively 
high, the azimuthal anisotropy also (v 2 ) starts increasing 
and therefore, increasing the background to CME. 


(b) 

FIG. 1. (a) (Color online) Body-tip configuration in the lab¬ 
oratory frame of reference. The impact parameter is along 
the X-axis, (b) The cross sectional view of the overlap region 
(shown by mesh) for a central (b = 0) body-tip collision. 

Due to the deformed shape, the U+U collisions can have 


unique orientations in which the magnetic field is very 
high in central collisions and as well as the azimuthal 
anisotropy is also very low. Therefore, U+U collisions 
may provide an unique opportunity to study these exotic 
effects in relativistic heavy-ion collisions. However, it has 
not been experimentally possible so far to unambiguously 
select specific configurations in U+U collisions. 

In this paper, we propose a methodology to select 
body-tip configuration from unbiased events in U+U col¬ 
lisions. The body-tip configuration is pictorially shown in 
Figure Oja), where the impact parameter ( b ) is along X- 
axis and the beam direction is along Z-axis. In this con¬ 
figuration, the (right going) uranium nuclei whose major 
axis is perpendicular to beam is called body and other 
one (left going) whose major axis is along the beam is 
called tip Q. As seen in Figure QJb), the overlap region 
in such a body-tip collisions is circular (shown by mesh). 
The nucleons which lies in the overlap are called partici¬ 
pants and those which lies outside the overlap region and 
do not take part in the collisions are called spectators. 
It is visible from Figure QJb) that one uranium nucleus 
get completely occluded into the other, leaving almost 
no spectators, where as the other one will always have 
some spectator from the non-overlapping regions. This 
gives rise to asymmetry in the spectator counts in the 
two opposite directions. We use this particular feature of 
this body-tip event configuration to separate it out from 
rest of all the other random configurations possible in 
deformed uranium nuclei. 


II. THE AMPT MODEL 

We have used A Multi Phase Transport, namely the 
AMPT model, for our analysis. U+U collisions are im¬ 
plemented in the AMPT model by deforming the Woods- 
Saxon profile Q as, 


1 + exp([r — R'}/a) 

R' = R[l + /3 2 Y 2 °(6)+faY°(0)\ 
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FIG. 2. (a) (Color online) Spectator neutron distribution for 0-5% central events in body-tip configuration of U+U collisions 
at snn = 193 GeV. (b) Energy deposited by spectator neutrons in ZDC for the same events of body-tip configuration as in 
(a), (c) Energy deposited by spectator neutrons in ZDC for all configurations in U+U collisions. The dotted lines show the 
selection range for body-tip events, (d) The spectator neutron energy in ZDC after selecting events using the dotted lines in 
(c). 


where p 0 is the normal nuclear density, R is the radius 
of the nucleus, and a denotes the surface diffuseness pa¬ 
rameter. Y) m ((9) denote spherical harmonics and 9 is the 
polar angle with respect to the symmetry axis of the 
nucleus. We have used the values of R, a , fa and f3i 
from [l]. The AMPT model, which is a hybrid transport 
model, has four main stages: the initial conditions, par- 
tonic interactions, the conversion from the partonic to the 
hadronic matter, and finally hadronic interactions It 
uses the same initial conditions from HIJING jf|. Scat- 
teringamong partons are modelled by Zhangs parton cas¬ 
cade Q , which calculates two-body parton scatterings us¬ 
ing cross sections from pQCD with screening masses. In 
the default AMPT model, partons are recombined with 
their parent strings and when they stop interacting, the 
resulting strings fragment into hadrons according to the 
Lund string fragmentation model 0- However in the 
string melting (SM) scenario, these strings are converted 
to soft partons and a quark coalescence model is used to 
combine parton into hadrons. The evolution dynamics of 
the hadronic matter is described by A Relativistic Trans¬ 
port (ART) model 0- We have used string melting 
(SM) mode of AMPT version 2.25t7, with parton-parton 
cross section of 10 mb which will give rise to substantial 
amount of Vi- A total of 0.6 million 0-5% central events 


were generated using this model for the analysis. 

III. RESULTS AND DISCUSSION 

In the experiments it is possible to get the measure of 
spectators using zero degree calorimeter (ZDC) detec¬ 
tor which lies very close to the beam pipe in the for¬ 
ward (and backward) direction. The ZDC detector gives 
an electrical signal which is proportional to the number 
of spectator neutrons. Therefore, we will use only the 
neutrons from the spectators for our study. In Figure [2] 
we show the spectator neutron correlation for both body 
and tip oriented Uranium nuclei in body-tip collisions. 
As seen from the figure, the spectator neutron counts 
are not symmetric for body-tip collisions. In a real ex¬ 
periment, these spectator neutrons are detected by Zero 
Degree Calorimeters (ZDC). Therefore, it is worthwhile 
to convert this spectator neutrons into practically mea¬ 
surable ZDC signal. We use the ZDC response of STAR 
experiment from Ref. pjj. From Fig.6 of Ref. 0] we 
find that the resolution of single neutron ZDC response 
is 18%. Therefore, we smear the energy deposited by each 
individual spectator neutron by a Gaussian distribution 
with the width of 18% of the mean value to imitate the 
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response in the ZDC detector. The mean value of energy 
deposited by a single neutron for this study is 96.5 GeV. 
The energy deposited by the spectator neutrons in the 
ZDC, event by event, is shown in Figure[5Jb). We can se¬ 
lect the body-tip events from all the other configurations 
using this correlation shown in Figure [2jb). One such se¬ 
lection procedure is shown by the dotted line, with slope 
= 0.25 and intercept = -180. We select all the events 
which lies below this line, therefore selecting the events 
with asymmetric spectator neutrons. The ZDC response 
for all possible configurations in U+U collisions is shown 
in Figure (He). Since both left going and right going nu¬ 
clei can be in either body or in tip orientation, therefore, 
we select these events along both (left and right) axis. 
Two dotted lines show the selection ranges for the possi¬ 
ble body-tip events. Figurc[2])d) shows the ZDC response 
of the selected events. 

One way to differentiate between the minimum bias (all 
possible configurations) and the body-tip configurations 
is to look at the variable S v which is defined as, 


S v 


YtVidN/dr]) 

Ntot 


(3) 


where, N to t is the number of particles within pseudora¬ 
pidity range, - 1.0 < 77 < 1.0 and the summation is over all 
particles in the event. Figure [3] shows the variable S v as a 
function of N to t for minimum bias, events selected with 
cuts on the ZDC signal from minimum bias configura¬ 
tion, and pure body-tip events. The S v for minmum bias 
configurations lies close to zero, suggesting symmetry in 
particle production. The particle production in body-tip 
events are asymmetric in 77 as shown by solid squares in 
Fig. [3] Whereas, when selecting body-tip events from 
minimum bias configurations, both the projectile and 
the target can either be in body or in tip configuration. 
Therefore we have two set of S v (when projectile is body, 
target is tip and vice-versa) for the selected events. The 
difference observed in S v of selected events and mini¬ 
mum bias configurations enhances the possibility of our 
method to select the body-tip events in real experiments. 

Now that we think we have potentially selected body- 
tip events from all configurations, we are going to look at 
V 2 of selected events. As the overlap region in a central 
body-tip collision is circular, therefore we expect that 
magnitude of the i >2 for selected events should be less 
compared to a set of events with an specific selection of 
configuration in U+U collisions. Figure 0 shows the V 2 
of charged particles in mid-rapidity (- 1 . 0 < 77 < 1 . 0 ), mea¬ 
sured with respect to participant plane (’Fpp) for 0-5% 
central events in U+U collisions at ^/snn = 193 GeV. 
The detailed method of calculating 'Fpp using the po¬ 
sition co-ordinates of participants (nucleons which take 
part in the collisions) is given in Ref. 0 ]. From Fig¬ 
ure UK a) we see that the V 2 in body-tip events is lower 
compared to minium bias configurations in U+U colli¬ 
sions. Now the selected body-tip events might also con¬ 
tain some small amount of other possible configurations 
of U+U collisions. Therefore, we have varied the slope of 


the dotted line shown in Figure 0c) to make the selec¬ 
tion more strict or relaxed. The corresponding change in 
the V 2 is also shown in the Figure 0(a). The magnitude 
of charged particle V 2 (pr) hi the selected body-tip events 
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FIG. 3. (Color online) S v (see Eqn. [3j) for minimum bias, 
selected body-tip events from minimum bias, and pure body- 
tip events. 

are systematically 25% lower than that in all configura¬ 
tions. The V 2 (pr) for purely body-tip events (i.e., orien¬ 
tation of Uranium nuclei are fixed according to body and 
tip orientation) in AMPT is also shown as solid curve in 
Figure HKa). The dotted and dashed lines in Figure 0 a) 
corresponds to the V 2 measured with respect to reaction 
plane angle (T r ) in body-tip and selected body-tip events 
from all configurations respectively. The difference in the 
UiK’I'r} and ^{pp} shown in Figure 0(a) arises due the 
fluctuations in the initial participant distribution. The 
reaction plane angle T r is defined as the angle between 
the impact parameter and the X-axis and is a known 
quantity in the AMPT model. In experiments, however, 
the position of the participant nucleons are unknown. 
In this scenario, the event plane is calculated using the 
anisotropic distribution of the produced particles 0 - 
We followed 77 sub-event plane method to calculate V 2 
of charged particles. In this method, each event is di¬ 
vided into two uncorrelated sub-events in two different 77 
windows. Then 2 nd order event plane (T 2 ) is calculated 
separately in both of these sub-events. Each particle is 
then correlated with the event plane of opposite 77 so as to 
remove the self-correlation effect 0 - The V 2 result ob¬ 
tained in this method is then corrected for the 77 sub-event 
plane resolution 0 ], We have followed event-by-event 
resolution correction jTt| for our analysis. The details of 
the procedure of event plane calculation and resolution 
correction can be found in 0 |. The U 2 {EP} results for 
0-5% central U+U collisions, calculated with event plane 
method, are shown in Figure 0b). The V 2 {EP} for all 
configuration is shown by open markers in Figure 0(b). 
The u 2 {EP} for different ZDC cuts are shown by solid 
markers. As seen in Figure0b), U 2 {EP} is also system- 
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atically 25% lower for selected body-tip events compared 
to u 2 {EP} of all configurations. The ^{EP} for purely 
body-tip events is shown by solid curve in Figure 0](b). 
As we see in Figure |4^b), with decreasing slope param¬ 
eter (i.e., selecting events with higher spectator neutron 
asymmetry), r; 2 {EP} of selected events systematically de¬ 
creases compared to that of all configurations. Therefore, 
using lower slope value (higher spectator asymmetry), 
we can enhance the selection of body-tip events from all 
configurations in U+U collisions. Figure 01(c) and Fig¬ 
ure 01(d) shows the distribution of the angle of the major 
axis of the Uranium with the beam axis for body (9 b ) and 
tip (9 t ) orientation for slope = 0.25. 9 b and 9 t for two 
different slope values and for all configuration in U+U 
collisions are also shown for comparison. Although lower 
slope values enhance the selection of body-tip events from 
all configuration, it also results in reduced event statis¬ 
tics. The selected events may possibly contain some other 
configurations. Therefore, we have calculated the pu¬ 
rity of the selected events. By purity we infer that how 
much of the selected events can be specified as body-tip 
type events. Therefore, purity(%) = (number of selected 
events with 9 t , 9 P within the limits of relaxation 9 R ) / (to¬ 
tal number of events satisfying the ZDC cut) x 100. Here 
we call an selected event to be body-tip if -9 R < 9 t < 9 r 
and (n/2 - 9 R ) < 9 b < (9 R + n/2). 


TABLE I. Purity (in %) of selected events for different slope 
parameter and angular relaxation (see text for details). 


Angular 
, relaxation 

Pe t ^\(0«) 
parameter 

±10° 

±20° 

±30° 

0.10 

34% 

62% 

74% 

0.15 

28% 

52% 

70% 

0.25 

26% 

48% 

67% 

0.35 

24% 

45% 

63% 


Since the angular distribution of both of the selected Ura¬ 
nium nuclei has some finite width, therefore, we assume 
that out of the selected events, the number of events for 
which 9 b and 9 t lies within zt9 R degree of the correspond¬ 
ing default values in body-tip events ( 9 b = tt/2 and 9 t = 


0) is the pure body-tip sample. It is worth mentioning 
that the purity of the selected events depends on the 
relaxation on the angular width (9 R ) we set to classify 
the selected event as body-tip. Hence we calculated the 
purity of the selected event sample for different angular 
relaxation and the results are listed in Table HI As seen 
from Table [U the purity of selected events increases as 
we decrease the value of slope parameter or if we increase 
the angular relaxation. As can be found in Tabic [TJ more 
than 70% purity can be achieved using low slope value 
(i.e., higher order of spectator neutron asymmetry). 

IV. SUMMARY 

We present an experimental procedure to select the body- 
tip configuration among all possible configuration in 0- 
5% central U+U collisions at yj snn = 193 GeV. We 
found that the spectator neutron energy deposited in the 
Zero Degree Calorimeter (ZDC) is a useful tool to select 
body-tip oriented events in central U+U collisions. We 
are able to select body-tip configuration with conditions 
applied on spectator neutron asymmetry simulated with 
the ZDC. We have used a new variable S v to differenti¬ 
ate between the body-tip and the minimum bias config¬ 
urations. Elliptic flow ( V 2 ) is calculated for the selected 
events with respect to both participant plane angle (d'pp) 
and event plane angle (\H 2 ). As expected, v 2 of selected 
events is found to be systematically lower compared to 
that in all configurations in U+U collisions. The ZDC 
selection cut (slope) was varied and it was found that se¬ 
lecting events with higher spectator neutron asymmetry 
results in lower v 2 values which tends to match with V 2 
of pure body-tip events. Finally we calculated the pu¬ 
rity of the selected events from all configurations in the 
U+U collisions. We observed that purity increases for 
decreasing slope parameter. In other words, if we ap¬ 
ply cut selecting higher order of left-right spectator neu¬ 
tron asymmetry, then the purity of the selected events 
increases. 
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FIG. 4. (Color online) (a) V2 {pp} for 0-5% central U+U collisions at v /snn = 193 GeV without ZDC cut (open markers) and 
with ZDC cut (solid markers), (b) V 2 {EP} for the same events without ZDC cut (open markers) and with ZDC cut (solid 
markers). Solid line in (a) and (b) corresponds to V2 PP and vep for pure body-tip events, (c) Distribution of the angle of the 
major axis with the beam axis for body (9b) and for (d) tip ( 9t ) oriented Uranium nuclei for different slope parameters. For 
pure body-tip events, (n/2 — 0.005) < 9b < (tt/ 2 + 0.005) and 0 < 9t < 0.05 (radian). 
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